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Abstract 

Let G = {V, E) be an undirected graph with a nonnegative edge- weight 
function w. The routing cost of a spanning tree T of G is dxiu, v), 

where dxi^u, v) denotes the weight of the simple u-v path in T. The Mini- 
mum Routing Cost Spanning Tree (MRCT) problem [WLB+00] asks 
for a spanning tree of G with the minimum routing cost. In this paper, we 
parallelize several previously proposed approximation algorithms for the 
MRCT problem and some of its variants. Let e > be an arbitrary con- 
stant. When the edge-weight function w is given in unary, we parallelize 
the (4/3 + e)-approximation algorithm for the MRCT problem [WCTOOb] 
by implementing it using an TZAfC circuit. There are other variants of the 
MRCT problem. In the Sum- Requirement Optimal Communication 
Spanning Tree (SROCT) problem [WCTOOa], each vertex u is associated 
with a requirement r{u) > 0. The objective is to find a spanning tree T of 
G minimizing ^,^y {r{u) + r{v)) driu, v). When the edge-weight func- 
tion w and the vertc^x-K^quirenient function r are given in unary, we paral- 
lelize the 2-approximation algorithm for the SROCT problem [WCTOOa] 
by realizing it using TZAfC circuits, with a slight degradation in the ap- 
proximation ratio from 2 to 2 + o(l). In the weighted 2- MRCT problem 
[Wu02], we have additional inputs si, S2 G V and A > L The objective is 
to find a spanning tree T of G minimizing J2vev ^^t{si,v) + dT{s2,v). 
When the edge-weight function w is given in unary, we parallelize the 
2-approximation algorithm [Wu02] into TZAfC circuits, with a slight degra- 
dation in the approximation ratio from 2 to 2 -|- o(l). To the best of our 
knowledge, our results are the first parallelized approximation algorithms 
for the MRCT problem and its variants. 

1 Introduction 

Let G = (y, E) be an undirected graph with a nonnegative edge-weight function 
w. The routing cost of a spanning tree T of G is driu, v) where dxiu, v) 

is the weight of any shortest u-v path in T, or cquivalently, the weight of the 
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simple u-v path in T. The Minimum Routing Cost Spanning Tree (MRCT) 
problem [WLB+00] asks for a spanning tree T of G with the minimum routing 
cost. It is also known as the Shortest Total Path Length Spanning Tree 
problem. The MRCT problem is first proposed by Hu [Hu74], who referred 
to the problem as the Optimum Distance Spanning Tree Problem. In 
Hu's formulation of the more general Optimum Communication Spanning 
Tree (OCT) problem [Hu74], an additional value Tu,v > is given for each 
pair {u,v) of vertices. The communication cost [Hu74] of a spanning tree T 
of G is Yluv€V ^t{u,v). The OCT problem asks for a spanning tree of G 
with the minimum communication cost. When G is a complete graph and the 
edge- weight function w obeys the triangle inequality, a randomized 0(log \ V\)- 
approximation algorithm is known for the OCT problem [Bar98, CCGG98, 
WLB+00, FRT03]. The MRCT problem is the special case of the OCT problem 
when Tu^v = 1 for all u,v £ V. 

The MRCT problem has applications in network design [Hu74, JLK78] as 
well as multiple sequences ahgnment in computational biology [FD87, Pev92, 
Gus93, BLP94, WLB+00]. Unfortunately, it is shown to be MV-hard [JLK78], 
and it is A^P-hard even when all edge weights arc equal [JLK78, GJ79] or when 
the edge- weight function obeys the triangle inequality [WLB+00]. 

Exact and approximation algorithms for the MRCT problem have been 
extensively researched [BFW73, Hoa73, DF79, Won80, WCTOOb, WLB+00, 
FLS02]. Boyce et al. [BFW73], Hoang [Hoa73] and Dionne and Florian [DF79] 
study branch- and-bound algorithms as well as heuristic approximation algo- 
rithms for the Optimal Network Design problem [BFW73], which includes 
the MRCT problem as a special case. Fischetti et al. [FLS02] give exact algo- 
rithms for the MRCT problem while avoiding exhaustive search. Wong [Won80] 
gives a polynomial-time 2-approximation algorithm for the MRCT problem. 
That is, he gives a polynomial-time algorithm that, given a graph G = {V, E) 
with a nonnegative edge- weight function w, outputs a spanning tree of G whose 
routing cost is at most 2 times the minimum. Subsequent work by Wu et al. 
[WCTOOb] shows a different polynomial-time 2-approximation algorithm as well 
as polynomial-time 15/8,3/2 and (4/3 -|- e)-approximation algorithms for the 
MRCT problem, where e > is an arbitrary constant. Their results are later 
improved by Wu et al. [WLB+00] to give a polynomial-time approximation 
scheme (PTAS) [CLRSOl] for the MRCT problem. That is, a polynomial-time 
(1 -|- e)-approximation algorithm is given for any constant e > 0. 

There are other variants of the MRCT problem that also have applica- 
tions in network design [WLB+00, WCTOOa, WCTOOc, Wu02]. In the SUM- 
Requirement Optimal Communication Spanning Tree (SROCT) prob- 
lem [WCTOOa], each vertex u is associated with a requirement r{u) > 0. The ob- 
jective is to find a spanning tree T of G minimizing ^^y{r{u)+r{v)) dxiu, v). 
The Product-Requirement Optimal Communication Spanning Tree 
(PROCT) [WCTOOa] problem is to find a spanning tree T of G minimizing 

vev ^(^) ''(^) '^Tiu, v)- The SROCT and PROCT problems are clearly gen- 
eralizations of the MRCT problem. 

Wu et al. [WCTOOa] give a 2-approximation algorithm for the SROCT 
problem. They also propose a 1.577-approximation algorithm [WCTOOa] for 



2 



the PROCT problem. The result is improved by Wu et al. [WCTOOc] to yield 
a polynomial-time approximation scheme (PTAS) for the PROCT problem. 

Another variant of the MRCT problem is the 2-MRCT problem [Wu02]. 
In this problem, except for G = {V,E) and w : E ^ Rq, we are given two 
source vertices si,S2 G V. The objective is to find a spanning tree T of G 
minimizing "^^^v ^^{''^1,^) + dT{s2,v). This problem is AAP-hard even when 
w obeys the triangle inequality [Wu02]. Wu [Wu02] shows a 2-approximation 
algorithm as well as a PTAS for this problem. A variant of the 2-MRCT 
problem is the weighted 2-MRCT problem [Wu02] where an additional A > 1 
is given as input. The objective is to find a spanning tree T of G minimizing 
X^„gy A (iT(si, u) + c?t(s2! ^^)- Wu [Wu02] proposes a 2-approximation algorithm 
for the weighted 2-MRCT problem. When the edge-weight function w obeys the 
triangle inequality, there is a PTAS for the weighted 2-MRCT problem [Wu02] . 

In this paper, however, we will focus on parallelizing the approximation 
algorithms for the above problems. We first describe our results concerning the 
MRCT problem. For an arbitrary e > and when the edge- weight function w 
is given in unary, we show that the (4/3-|-e)-approximation algorithm proposed 
by Wu ct al. [WCTOOb] can be implemented by an TZAfC circuit. That is, the 
approximation algorithm can be performed by a uniform polynomial-size circuit 
[Pap94] with random gates and poly-logarithmic depth. Indeed, with a small 
probability our parallelized algorithm may fail to find a (4/3-1- e)-approximate 
solution, in which case it outputs "fail." Thus, our algorithm does not fail (to 
find a (4/3-|-e)-approximate solution) without ever knowing that it fails, which 
is a desirable property for randomized algorithms with a small probability of 
failure. 

We now turn to describe our results concerning the SROCT problem. When 
the edge-weight and the vertex-requirement functions are given in unary, we 
parallelize the 2-approximation algorithm [WCTOOa] by realizing it using TZAfC 
circuits, with a slight degradation in the approximation ratio (from the currently 
best 2 to our 2 -|- o(l)). Still, with a small probability our algorithm may fail 
to output a (2 -f o(l))-approximate solution, in which case it knows the failure 
and outputs "fail." 

Finally, for the weighted 2-MRCT problem with the edge-weight function 

given in unary, we parallelize the 2-approximation algorithm [Wu02] into TZJ\fC 
circuits, with a slight degradation in the approximation ratio (from the currently 
best 2 to our 2 -|- o(l)). Again, there is a small probability that our algorithm 
fails to find a (2 -|- o(l))-approximate solution, in which case it outputs "fail." 

To the best of our knowledge, our results are the first efforts towards paral- 
lelized approximation algorithms for the MRCT problem and its variants. Our 
results open up new opportunities to compute approximate solutions to the 
above problems in parallel poly-logarithmic time. In the applications of the 
MRCT problem to network design [Hu74, JLK78] as well as the applications 
of the SROCT and PROCT problems to network design [WCTOOa, WCTOOc], 
the network is often modeled as a graph with a nonnegative edge-weight func- 
tion representing the distances between pairs of nodes. Although approximate 
solutions to the aforementioned problems (MRCT, SROCT, PROCT, weighted 
2-MRCT) are attainable in polynomial time, in any real networking environ- 
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ment, however, the cost of traffic between any pair of nodes may vary over 
time. Tlius, it is highly desirable to be able to compute approximate solutions 
to these problems as fast as possible, so as to reduce the risk that the traffic 
costs change during the computation. Our results imply that approximate so- 
lutions to the MRCT, SROCT and weighted 2-MRCT problems can indeed by 
computed in parallel poly-logarithmic time on multiprocessors. 

For other applications of the MRCT problem where the data does not change 
quickly over time, for example multiple sequences alignment in computational 
biology [FD87, Pev92, Gus93, BLP94, WLB+00], being able to compute ap- 
proximate solutions to the MRCT problem in parallel sublinear time is still 
beneficial. Indeed, Fischer [FisOl] argues that in many practical applications 
today, the input size is so large that even performing linear-time computations 
is too time-consuming. Certainly, multiple sequences alignment in computa- 
tional biology constitutes a good example where the input size is usually too 
large. It is therefore a desirable property that our algorithms operate in parallel 
sublinear time, and in fact poly- logarithmic time. 

The main idea underlying our proofs is that many of the previously proposed 
approximation algorithms for the MRCT, SROCT and weighted 2-MRCT prob- 
lems rely heavily on finding shortest paths between pairs of vertices in a graph. 
This motivates applying the well-known result that MC C J\fC to parallelize 
these algorithms since we can guess a path (possibly the shortest one) between 
two vertices of a graph in nondeterministic logarithmic space. There is the com- 
plication that, in our proofs, we will often need to generate the same shortest 
path between two vertices u, v, whenever a shortest u-v path is needed. For 
this purpose, we use the isolation technique [Wig94, GW96, RAOO] to slightly 
modify the edge-weight function of the input graph, so that there is exactly 
one shortest path between each pair of vertices with high probability. We then 
apply the double-counting technique [RAOO] to decide whether the input graph 
(with the modified edge-weight function) exhibits a unique shortest path be- 
tween each pair u, v of vertices. If so, we are able use the double counting 
technique to generate the unique shortest u-v path whenever it is needed. The 
whole procedure runs in unambiguous logarithmic space and our results follow 
by UjC C MjC C J\fC. The approximation ratio would be slightly degraded. The 
degradation comes from randomly modifying the edge- weight function when we 
apply the isolation technique. 

Our paper is organized as follows. Section 2 provides the basic defini- 
tions. Section 3 presents the parallelized (4/3 -|- e)-approximation algorithm 
for the MRCT problem. Section 4-5 describe our parallelized approximation 
algorithms for the SROCT and the weighted 2-MRCT problems, respectively. 
Section 6 concludes the paper. Proofs are given in the appendix for references. 

2 Notations and basic facts 

Throughout this paper, graphs are simple undirected graphs [WesOl]. That is, 
we disallow parallel edges and self-loops. There will always be a nonnegative 
edge-weight function mapping each edge to a nonnegative real number. For a 
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graph G, V{G) is its vertex set and E{G) is its edge set. Let ii be a subgraph 
of G. A path P connects a vertex v to R (or V{R)) if one endpoint of P is f 
and the other is in V{R). An edge connecting two vertices u and v is denoted 
uv. A path connecting two vertices u and v is said to be a u-v path. A path 
{vq, . . . ,Vk) is one which traverses vq, ■ ■ ■ ,Vk, in that order. A simple path is 
a path that traverses each vertex at most once [WesOl]. A graph G contains 
another graph G' if G' is a subgraph of G. The set of nonnegative real numbers 
is denoted Mq". 

Definition 1. Let G = {V,E) be an undirected graph and w : E ^ Mq be 
a nonnegative edge-weight function. The lexicographical ordering on V is that 
of the encodings of the vertices in V, assuming any reasonable encoding of a 
graph. Let u,v £ V and R be a subgraph of G. The sum of edge weights of 
R is denoted w{R). When R is a path, w{R) is called the weight or length of 
R. For x,y G V, we use dcix, y) to denote the weight of any shortest x-y path. 
We use daix, R) ( or doix, V{R)) ) for min2gv'(^) dc^x., z). The lexicographically 
first vertex x' G V{R) satisfying dG{x,x') = dG{x,R) is denoted closest{x, R) 
(or closest{x, V{R)) ). The set of all shortest paths connecting u and v is denoted 
SPg{u,v). SPg{u,R) (or SPg{u,V{R))) is the set of shortest paths connecting 
u and R. That is, SPg{u,R) is the set of paths that connect u and R and have 
weight equal to dG{u,R). For k £N, Sk^u denotes the set of vertices reachable 

(k) 

from u with at most k edges. SPu^i, denotes the set of all shortest paths among 

(k) 

those u-v paths with at most k edges. That is, SPu^l is the set of u-v paths with 
at most k edges whose weight is not larger than any other path with at most k 
edges. The union of two graphs Gi = {yi,E{) and G2 = {V2,E2) is the graph 
(Vi UV2,Ei U E2). The graph G is strongly min-unique with respect to w if for 



all k and u,v ^ V, we have 



< 1. When w is clear from the context. 



we may simply say that G is strongly min-unique without referring to w. 



if 



In Definition 1, it is not hard to show that G = (V, E) is strongly min-unique 
SpI^I < 1 for all A; G {0, ... , \V\ - 1} and u,v eV, provided \V\ > 3. 
The MRCT of a graph, standing for its Minimum Routing Cost spanning 
Tree, is defined below. 

Definition 2. ([WCO4]) Given a connected graph G = (y,E) with a nonneg- 
ative edge-weight function w : E ^ Mq , the routing cost Cw{T) of a spanning 
tree T of G is '^I'iu^v) . A spanning tree of G with the minimum routing 

cost is an MRCT of G, which is denoted by MRCT{G) for convenience. 

The MRCT problem asks for MRCT(G) on input G, w. The following fact 
shows that the routing cost of a tree can be computed efficiently. 

Fact 3. ([WCO4]) Let G be a graph with a nonnegative edge-weight function 
w : E{T) M.Q and T be a spanning tree of G. For each edge e G E{T), let 
Te^i and Te^2 be the two trees formed by removing e from T. We have Cw(T) = 
E'eeE(T) 2|V^(Te,i)| \V {T,,2)\ w{e) and c^T) < \V{T)\y2 ■ nm^.^EiT) w{e). 

To ease the description, we introduce the following definition. 
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Definition 4. ([RAOO]) A nondeterministic Turing machine M outputs a string 

s unambiguously on input x if it outputs s on exactly one non-rejecting com- 
putation branch, and rejects x on all other computation branches. The unam- 
biguously output string s is also denoted M{x). 

Throughout this paper, when a nondeterministic Turing machine A runs 
or simulates another nondeterministic machine B, it means that A runs B 
and make nondeterministic branches as B does. It does not mean that A 
enumerates all computation branches of B and simulate them deterministically. 
For convenience, A does not necessarily have to output B^s output. Instead, it 
may extract portions of B's output for output. 

We will need the notion of a general star to introduce the approximation 
algorithms for the MRCT problem. 

Definition 5. ([WCO4]) Let G be a connected graph with a nonnegative edge- 
weight function w : E ^ and S be a subtree of G. A spanning tree T 
containing S is a general star with core S if each vertex u E V satisfies 
dT{u,S) = dc{u,S). When V{S) = {v} is a singleton, a general star with 
core S is also called a shortest path tree rooted at v. 

Given any subtree S oi G = {V,E), a general star with core S exists 
[WCTOOb]. This follows by observing that for any shortest path P connect- 
ing u € V and S, the part of P from any vertex x G V{P) to S constitutes a 
shortest path connecting x and S. 

The notion of a metric graph is defined below. 

Definition 6. ([WLB^ 00]) A complete graph G with a non-negative edge- 
weight function w is metric if w{xy) + w{yz) > w{xz) for all x,y,z & V{G). 

3 A parallelized (4/3 + e) -approximation for MRCT 

We begin with the following form of the famous isolation lemma. It is implicit 
in some previous works [Wig94, GW96, RAOO]. 

Theorem 7. ([RAOO]) Let G = {V,E) be a graph with a nonnegative edge- 
weight w : E ^ Mq". Let Wr ■ E ^ Mq" assign the weight of each e & E 
independently and randomly from the uniform distribution over a set W C M+ 
With probability at least 1 — |y|^/(2|VF|), the graph G is strongly min-unique 
with respect to w -\- Wr- 

The following theorem is implicit in [RAOO]. It uses the double counting 
technique [RAOO] similar to the inductive counting technique used to prove the 
Immerman-Szelepcsenyi theorem [Imm88, Sze88]. 

Theorem 8. ([RAOO]) There is a nondeterministic logarithmic- space Turing 
machine FIND-PATH that, on input a graph G = {V,E) with a nonnegative 
edge-weight function w : E ^ {0, . . . , poly{\V\)} and two vertices s,t G V, 
satisfies the following conditions. 
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1. If G is not strongly min-unique, then FIND-PATH outputs "not strongly 
min-unique" unambiguously. 

2. IfG is strongly min-unique and has an s-t path, then FIND-PATH outputs 
the unique path P G SPG{s,t) and its weight w{P) unambiguously. The 
edges in P are output in the direction going from s to t. 

3. If G is strongly min-unique and does not have an s-t path, then FIND- 
PATH has no accepting computation branches. 

The following theorem is due to Wu et al. [WCTOOb] . 

Theorem 9. ([WCTOOb]) Let r e N be a constant and G = {V,E) be a 

connected, strongly min-unique graph with a nonnegative edge-weight function 
w : E ^ R+. For 1 < k < r + 4 and S = (vi, ...,Vk) G V'', let Ri^s be the 
subgraph of G containing only vi. For 2 < i < k, let Ri^s = Ri-i,s U Pi^s where 
Pi^S € SPq {vi, closest{vi, Ri-i^s))) is the unique shortest path connecting Vi and 
closest{vi, Ri-i^s)- For some 1 < k <r-\-A and S eV'', every general star T 
with core Rk,s satisfies 

c^{T) <{^- + c^{MRCT{G)). 

That Ri^s in Theorem 9 is a subtree of G for 2 < z < A; is easily shown 
because w{e) > for each e £ E hy the strong min-uniqueness of G. The 
core Rk,s in Theorem 9 is unambiguously computable in logarithmic space on 
strongly min-unique connected graphs. To show this, we need the following 
lemma. 

Lemma 10. There is a nondeterministic logarithmic-space Turing machine 
ADD-PATH that, on input a strongly min-unique connected graph G = (V, E) 
with a nonnegative edge-weight function w : E ^ {0, . . . , po/j/(|F|)}, a subgraph 
R of G and a vertex v G V, outputs the unique path P G SPq {v, closest{v, R)) 
unambiguously. 

With Lemma 10, we are able to compute the core Rk,s in Theorem 9 un- 
ambiguously in logarithmic space on strongly min-unique connected graphs. 

Lemma 11. Let r gN be a constant. There is a nondeterministic logarithmic- 
space Turing machine CORE that, on input a strongly min-unique connected 
graph G = {V, E) with a nonnegative edge- weight function w : E ^ {0, . . . , poly{\V\)} 
and S = {vi, . . . ,Vk) G where 1 < < r + 4, unambiguously outputs Rk,s 
defined below. Ri^s is the subgraph of G containing only vi. For 2 < i < k, 
Ri,S = Ri-i,S U Pi,s where 

Pi^s e SPg {vi, closest{vi, Ri-i,s)) 

is the unique shortest path connecting vi and closest{vi, Ri-i^s)- 

With Theorem 9 and Lemma 11, it is not hard to show the following fact. 
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Fact 12. Let r £ N be a constant and G = {V, E) be a strongly min-unique, con- 
nected graph with a nonnegative edge-weight function w : E ^ {0, . . . ,poly{\V\)}. 
For a sequence S of at most r + 4 vertices in V, let Cs = CORE{G, w, S) and 
Pu e SPg {u, closest{u, Cs)) forueV\ V{Cs). Then 

Ts = CsU U ^« 

u&V\V{Cs) 

is a general star with core S, and 

cUTs) <(^l + ^^4^) • cUMRCT{G)) 

for some S. 

The general star with a core in Fact 12 can be computed unambiguously in 
logarithmic space, as the next lemma shows. 

Lemma 13. Let r £ N be a constant. There is a nondeterministic logarithmic- 
space Turing machine STAR that, on input a strongly min-unique connected 
graph G = {V, E) with a nonnegative edge-weight function w : E ^ {0, . . . , poly{\V\)} 
and a sequence S of at most r + 4 vertices in V, outputs Cs = CORE{G, w, S) 
and each unique path in SPq {u, closest{u, Cs)) for u E V \ V{Cs) unambigu- 
ously. 

The following lemma allows unambiguous logarithmic-space computation of 
the routing cost of a tree. 

Lemma 14. There is a nondeterministic logarithmic-space Turing machine 
ROUT-PAIR that, on input a tree T with a nonnegative edge-weight function 
w : E{T) — >■ {0, . . . , poly{\V{T)\)} and s,t G V{T), unambiguously outputs the 
unique simple path P* connecting s and t in T and w{P*). 

Combining Fact 12 and Lemmas 13-14 gives the following lemma. 

Lemma 15. Letr gN be a constant. There exists a nondeterministic logarithmic- 
space Turing machine APPROX that, on input a strongly min-unique con- 
nected graph G = (F, E) with a nonnegative edge-weight function w : E ^ 
{0, . . . ,poZy(|F|)}, unambiguously outputs a spanning tree T of G with 

c^T) <(^^ + • c^{MRCT{G)). 

The following lemma will be useful. 

Lemma 16. Let a > be a constant. Let G = {V,E) be a graph with a 
nonnegative edge-weight function w : E ^ Wq, and the minimum nonzero 

weight assigned by w, if it exists, is at least 1. Let Ti and T2 be spanning trees 
of G. Let Wr assign to each edge e E E a nonnegative weight w{e) < 1/\V\'^ and 
w' = w + Wr. Then 

Cyj'iTi) <ac^'{T2) (1) 
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implies 



cUTi)<a{l + 



1 



) 



(2) 



2\V\ 



for sufficiently large \ V\. 

Combining Theorem 7 and Lemma 15-16 yields the following theorem. 

Theorem 17. Let e > be a constant. There is an nMC'^ algorithm PARAL- 
LEL that, on input a weighted undirected graph G = {V, E) with a nonnegative 
edge-weight function w : E ^ {0, . . . , poly{\V\)}, satisfies the following. 

1. If G is disconnected, then PAR ALLEL{G,w) outputs "disconnected." 

2. If G is connected, then PARALLEL(G , w) outputs a spanning tree of 
G unambiguously or outputs "fail" unambiguously. The probability that 
PARALLEL{G) outputs a spanning tree of G unambiguously is at least 
1 - 1/(2|F|). If PARALLEL{G) outputs a spanning tree T of G unam- 
biguously, then 



We begin this section with the following definition. 

Definition 18. ([WCTOOa, Wu02]) Let G = {V,E) be a graph with a non- 
negative edge-weight function w : E ^ and r : V ^ M.q be a requirement 
function on vertices. Let si,S2 G V be two vertices of G and T be a spanning 
tree of G. The sum-requirement communication (s.r.c.) cost ofT is 



The Sum-Requirement Optimal Communication Spanning Tree (SROCT) 
problem is to find a spanning tree T of G with the minimum value of Cw\t) 
over all spanning trees of G. We use SROCT{G) to denote an arbitrary span- 
ning tree of G with the minimum s.r.c. cost. The two-source routing cost ofT 
with sources si,S2 is 



The 2-MRCT problem is to find a spanning tree TofG with the minimum value 

(2) 

of Civ {T -I 81,82) over all spanning trees of G (in this problem si and S2 are part 
of the input). We use 2-MRCT{G) to denote an arbitrary spanning tree of G 
with the minimum two-source routing cost when the sources s\ , S2 are clear from 




4 The SROCT problem 



c'i\T)= J2 {r{u) + r{v))dT{u,v). 



c'i^HT,si,S2) = {dT{si,v) + dT{s2,v)) . 



9 



The weighted two-source routing cost of T with sources 

c!S\T,si,S2,X) 

= {XdT{si,v) + dT{s2,v)) . 
vev 

The weighted 2-MRCT problem is to find a spanning tree T of G with the mini- 
mum value of Cw\t, si, S2, A) over all spanning trees of G (in this problem si, S2 
and X are part of the input). We use W-2-MRCT{G) to denote an arbitrary 
spanning tree of G with the minimum weighted two-source routing cost when X 
and the sources si,S2 are clear from the context. 

The SROCT, 2-MRCT and weighted 2-MRCT problems are all A/':P-hard, 
even on metric graphs [WLB+00, WCTOOa, WCTOOc, Wu02]. 

The following theorem gives a 2-approximation solution to the SROCT prob- 
lem. 

Theorem 19. ([WCTOOa]) Let G = iV^E) be a connected graph with a non- 
negative edge-weight function w and a nonnegative vertex-requirement function 
r. There exists a vertex x E V such that any shortest path tree T rooted at x 
satisfies 

4^)(r) < 2cl^\SR0CT{G)). 

Theorems 7-8, 19 and Lemma 14 give the following parallelized 2-approximation 
solution to the SROCT problem. 

Theorem 20. There is an TZAfC"^ algorithm PARALLEL-SROCT that, on in- 
put a connected graph G = (V, E) with a nonnegative edge-weight function 
w : E {0, . . . ,poly{\V\)} and a nonnegative vertex-requirement function 
r : V ^ {0, . . . ,poZy(|F|)}, outputs a spanning T of G with 

c^^\T)<{2 + o{l))cl^\SROCT{G)) 

with high probability. If PARALLEL-SROCT does not output such a spanning 
tree, it outputs "fail. " 

5 Weighted 2-MRCT problem 

For the weighted 2-MRCT problem, we can assume without loss of generality 
that the two sources si,S2 are such that £^^(•81,52) > 0, where G is the input 
graph. Otherwise, the problem reduces to finding a shortest path tree rooted 
at si, which was implicitly done in the proof of Theorem 20. Wu [Wu02] has 
the following 2-approximation solution for the weighted 2-MRCT problem. 

Theorem 21. ([Wu02]) Let G = {V,E) be a connected graph with a nonnega- 
tive edge-weight function w : E —>■ M^, two sources si,S2 G V with dG{si,S2) > 
and A > 1. Denote 

Di{v) = (A -I- 1) dciv, si) -\- daisi, S2) 



the context. Let A > 1. 
si,S2 and weight X is 
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and 

D2{v) = (A + 1) daiv, S2) + A dcisi, §2) 

for veV. Let Zf = {v \ Di{v) < D2{v)} and Zlf = V\Z^. Let Q e SPg{si,S2) 
be arbitrary. Denote 

Q = {qo = si, . . . , qj,qj+i, ...,82) 

where qj+i is the first vertex on Q (in the direction from si to S2) that is not 
in Zf (it is easy to see that si G Zf). For each v e V, let Py^si G SPg{v,si) 
and Py^s2 e SPg{v, S2) be arbitrary. If Ti = Ut,ezj" P'o,si ^2 = [jy^z^ Pv,s2 
are trees, then T = T1UT2U qjqj+i is a spanning tree of G and 

c^^\T) < 2c^^\W-2-MRCT{G)). 

Theorems 7-8 and 21 and Lemma 14 yield the fohowing theorem. 

Theorem 22. There is an TZMC^ algorithm WEIGHTED-2-MRCT that, on 
input a graph G = {V,E) with a nonnegative edge-weight function w : E ^ 
{0, . . . ,poZ2/(|F|)}, s\,S2 &V and A > 1, with high probability outputs a spanning 
tree T with 

42)(r) < {2 + o{l))c'^;^\W-2-MRCT{G)). 
If WEIGHTED-2-MRCT does not output such a spanning tree, it outputs "fail. " 

We make the following concluding remark. All our algorithms are shown to 
be 7?,AAC^-computable by showing that they run in unambiguous logarithmic 
space and succeed in giving an approximate solution when the random input 
specifies an edge-weight function Wr such that G is strongly min-unique with 
respect to u; + ty^- By a method similar to that in [RAOO], we can also turn the 
random weight assignment into polynomially long advices. This is summarized 
below. 

Corollary 23. Let e > be a constant. There are II C/ poly algorithms for 
{4:/ 3+e) -approximating the MRCT problem, {2+ o{l)) -approximating the SROCT 
problem and {2 -\-o{l))- approximating the weighted 2-MRCT problem, where the 
respective edge-weight and vertex-requirement functions are given in unary. 

6 Conclusion 

We have given parallelized approximation algorithms for the minimum routing 
cost spanning tree problem and some of its variants. Our results show that, by 
exhibiting multiple processors, we can compute approximate solutions to the 
considered problems in parallel poly-logarithmic time. We hope this will shed 
hght on the many areas in which the considered problems are concerned, for 
example network design [Hu74, JLK78] and multiple sequences alignment in 
computational biology [FD87, Pev92, Gus93, BLP94, WLB+00]. 
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Appendix 

Proof of Theorem 7. The theorem is clearly true for \V\ < 2. If G is not strongly 
min-unique with respect to w + Wr and \V\ > 3, we have seen that there exist 



> 2 where the edge weights are 

le existence of an edge e € E such 



< k < \V\ - l,s,t e V such that 
given with respect to w + Wr- This implies t' 

(k) 

that at least one path in SP^ / contains e, and at least one does not. In this 
case we say that {k, s, t) blames e. Thus, the probability that G is not strongly 
min-unique is at most the sum over < k < \ V\ — 1, s,t e V and e G £^ of the 
probability that {k, s, t) blames e. 

For any k G {0, . . . , |V"| — 1}, s,t £ V,e E E and any partial weight assign- 
ment of Wr to E \ {e}, there is at most one assignment of Wr to e to make 
{k,s,t) blame e. This is because if {k,s,t) blames e when Wr{e) is assigned a 
value Wr,ej then increasing or decreasing the value of Wr{e) forces all shortest 
paths among those s-t paths with at most k edges to exclude or include e, 
respectively, making e no longer blamed by {k,s,t). Therefore, 

^Pr [(A;, s,t) blames e] 

0<k<\V\-ls,teV eeE 

\w 



< 



\W\ 

0<k<\V\-l s,teV eeE ' ' 

2\wy 



completing the proof. □ 
Proof of Theorem 8. Before describing how FIND-PATH works, we describe a 

(k) (k) 

few procedures that are useful. For k eN, let Ck = \Sks\- Let PrJ G SP^Cj be 
arbitrary and = ^^^^^ w{Ps'^v)- Note that the definition of does not 

(k) (k) 

depend on exactly which path in SP^,^ is chosen as Ps,v ■ It is clear that cq = 1 
and Eq = 0. 

We first introduce a nondeterministic logarithmic-space subroutine OUT- 

(k) 

PUT that outputs Sk,s unambiguously, given G, w, s, Cfe, and that SPs^t, = 1 
for each v G S^^s- OUTPUT just needs to nondeterministically guess each ver- 
tex X to be in or out of Sk^s, and if the guess is x G Sk.s then it outputs x. It 
verifies each guess of x G S^^g by nondeterministically guessing an s-x path with 
at most k edges and rejecting if it fails. Along the way OUTPUT counts the 
number c'j, of vertices verified to be in Sk,s and accumulates the weights of the 
guessed s-x paths (for x verified to be in Sk^s) in a variable S^. It then rejects if 
c'k / Cfc or E'^ / Efc. Clearly, guessing any vertex out of Sk,s to be in Sk,s results 
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in rejection. For a computation branch of OUTPUT not to reject, it must have 

reach c^., which requires successfuUy guessing an s-x path with at most k 
edges for eacli x G Sk.s- But to have S'^ not exceed Sfc, the guessed s-x path for 

each X £ Sk^s should be the unique one in SPi|2- So OUTPUT(G, ti;, s, c^, S^) 
has a unique non-rejecting computation branch, on which it correctly guesses 
whether each vertex x belongs to S^^g £^nd if so, correctly guesses the unique 
path in SpS. 



We now describe a procedure INDUCTIVE that computes Ck+i and S 



fc+i, 



and determines whether 
Ck and Sfc and that 



S'O 



> 1 for some x €z V unambiguously, given 

1 for each v £ Sk^s- For each vertex x G V, 
INDUCTIVE runs OUTPUT(G, tt;, s, Cfe, S^) to determine whether x G Sk.s 

(k) 

unambiguously and if so, accumulates the weight Wx of the unique path in SPs,:e 
as it is guessed by OUTPUT. For each u such that ux is an edge, INDUCTIVE 
also determines whether u G Sk,s unambiguously and if so, accumulates the 

(k) 

weight Wu of the unique path in SP^,^ as it is guessed. INDUCTIVE then 
computes the weight of any shortest s-x path with at most k + 1 edges as 



minQ 



where 



if X G -Sfe s and 



Q = {wx} U {wu + w{ux) I u G Sk,s, ux G E} 



Q = {wu + w{ux) I u G Sk,s, ux G -E} 

otherwise. If Q = 0, INDUCTIVE knows that x ^ Sk+\,s- Otherwise x G Sk+i,s 
and the weight of any path in SPif^^^ is known to be minQ. In case of a tie 



when computing minQ, INDUCTIVE knows that 
above for all x G F allows INDUCTIVE to compute Ck+i 



s.x 



> 1. Doing the 
'S'a:+i,s| and Sfc+i 



and determine whether 



op(fc+l) 
iJ-'^ s,x 



> 1 for some x V unambiguously. During 
the computation of minQ, INDUCTIVE does not store the set Q. Instead, 
INDUCTIVE computes the elements of Q one by one and stores the smallest 
element in Q that has been computed so far, as well as a flag indicating whether 
min Q is achieved by two elements at any time. 

We are now ready to describe how FIND-PATH works. Assume \V\ > 3. 
FIND-PATH starts with ccSq and repeatedly simulates INDUCTIVE un- 



Qp(^) 
'JJ^ S,X 



> 1 for some 



til it computes C|y|,E|y| or until it determines that 

< A; < |F| — 1 and x e V. Doing the above with each other vertex s' e V 
replacing the role of s guarantees that if G is not strongly min-unique, then 
we must find SP^!^|j. > 1 for some < k < \V\ ~ 1 and s',x G V. Instead, 
if G is strongly min-unique then FIND-PATH will compute all the way from 
Co, So to C|y|,i;|y|. It then runs OUTPUT(G, s, C|y|_i, E|y|_i). As we have 
seen, OUTPUT(G, w, s, c\v\-i,'^\v\-i) lias a unique non-rejecting computation 
branch, on which the unique shortest s-t path P G SP^'^'~^^ is correctly guessed 



13 



by OUTPUT. The weight w{P) is accumulated along the way. The strong min- 
uniqueness of G guarantees that P is also the unique path in SPg(s, t). □ 

Proof of Lemma 10. For each x G V{R), ADD-PATH runs FIND-PATH(G, w, x, 
to unambiguously generate the unique shortest path Px.v G SPG'(a;,f) and its 
weight w{Px,v)- In this way, ADD-PATH could compute minj.gy(/j) w[Px,v) as 
well as closes't(i;, R). Then ADD-PATH simulates FIND-PATH (G, w, v, closest(i;, 
to output the unique path P in SPg {v, closest('y, R)) unambiguously. □ 

Proof of Lemma 11. Clearly, CORE could output Ri^s unambiguously. Let 
2 < i < A;. To output Rj^s = Rj-i,S U Pj,s unambiguously, CORE recursively 
outputs Rj^i^s and then runs ADD-PATH(G, it;, i?j_i^5, f j) to unambiguously 
output Pj,s- There is an additional complication that CORE does not store 
Rj-i,S before calHng ADD-PATH. Instead, whenever ADD-PATH wants to read 
any bit encoding Rj^i^s, CORE recursively outputs Rj-i,s unambiguously on 
the fly to support the required bit. Each level of recursion uses up logarithmic 
space and the depth of recursion is at most r + 4, a constant. The space 
requirement is therefore logarithmic. 

□ 

Proof of Lemma 13. STAR begins by running CORE(G, w, S) to output Cs 
unambiguously. For any -u G V, STAR runs CORE(G, 5) to unambigu- 
ously determine whether u G V{Cs)- If u ^ V{Cs), STAR needs to out- 
put the unique path in SPg («, closest('U, C5)) . For this purpose, it computes 
closest(u, C5) as follows. For each w G V, STAR tests if w G V{Cs), again by 
running CORE(G, w;, S). If w ^ V{Cs), STAR goes on with the next v e V. 
Otherwise, STAR invokes FIND-PATH(G, ti;, i;) to generate the unique path 
Pu,v £ ^^g{u,v) and its weight w{Pu,v) unambiguously. STAR records the 

V G V{Cs) that has generated the smallest value of w{Pu,v) so far, favoring 
lexicographically smaller values of v in case of a tie. In the end, the recorded 

V must be closest ( n, C5) by the definition of closest(u, C5). At this time 
STAR just invokes FIND-PATH (G, ly, closcst(n, C5)) to output the unique 
path in SPg («, closest(ii, C^)) unambiguously. Doing the above for all -u G F 
does the job. □ 

Proof of Lemma I4. li s = t the task is trivial. We assume otherwise. ROUT- 
PAIR nondeterministically guesses a path P that does not enter a vertex im- 
mediately after it has left that vertex. If P is an s-t path, then ROUT-PAIR 
accepts, otherwise it rejects. The simple s-t path in T is the only s-t path that 
does not enter a vertex immediately after leaving it. Its weight w{P*) can be 
accumulated as it is guessed. □ 

Proof of Lemma 15. For each sequence S of at most r-|-4 vertices in V, Lemma 13 
enables us to unambiguously output Cs = CORE(G, w, S) and then each unique 
path Pu in SFq {u,closest{u,Cs)) for n G V\ V{Cs). Furthermore, Fact 12 
guarantees that 

Ts = CsU U Pu 

uev\v{Cs) 
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is a spanning tree of G and satisfies 

c^(Ts) < Q + ■ c^(MRCT(G)) 

for some S. Thus, we need only compute Cyj{Ts) unambiguously for each se- 
quence S of at most r + 4 vertices, and output Tg* unambiguously for the se- 
quence S* of at most r+4 vertices satisfying Cw{Ts*) = laing^yk^Kk^r+A Cw{Ts)- 
By Definition 2, Cw(Ts) can be computed unambiguously by running R0UT-PAIR(T5, w, s, t) 
for all pairs s,t G V and summing up the weight of the simple paths as they 
are output. There is a complication that APPROX does not store T5 before 
calling ROUT-PAIR. Instead, when ROUT-PAIR wants to read any bit in the 
encoding of Ts, APPROX runs STAR(G, w, S) to generate the required bit un- 
ambiguously on the fly. This enables us to obtain S* unambiguously and thus 
Ts* unambiguously by running STAR(G,i(;, 5*). □ 

Proof of Lemma 16. It is clear that either Cyj{T2) = or Cy,{T2) > 1. Also, 
Fact 3 implies that 

c..(T,) < ^. (3) 

If Cw{T2) = 0, Eq. (3) implies that 

Cw'{T2) = c^,(r2) < ^ 

and thus Cyji(Ti) < 1 for sufficiently large \V\ by Eq. (1). This implies Cw{Ti) < 
1 and thus Cw{Ti) = 0, establishing Eq. (2). 
If Cy,{T2) > 1, then Eq. (1) and (3) imply 

(Ti) < Cw'{Ti) 

< «(i + ^)c«,(r2). 

□ 

Proof of Theorem 17. We will show that PARALLEL needs only take a poly(|l^|)- 
long random input and do the rest of the computation in unambiguous logarith- 
mic space. The standard proof technique for showing that UC C J\f£. C J\fC'^ 
[Pap94, Sip05] then completes the proof. 

PARALLEL tests the connectedness of G by testing each pair of vertices 
for connectedness in logarithmic space [Rei05]. 

Below we assume that G is connected. If we assume that this theorem 
is true when w is not identically zero, then PARALLEL can also deal with 
the identically zero case by using the unit edge-weight function instead. The 
output spanning tree would have zero routing cost under the identically zero 
edge-weight function, so item 2 is still satisfied. Thus, we can assume without 
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loss of generality that w is not identically zero. Furthermore, we can normalize 
w to give mm^^E,w(e)7tow{e) > I- 

The random input to PARALLEL determines an edge-weight function Wr '■ 
E — > Mq" where for each e e E, Wr{e) is independently and randomly chosen 
from the uniform distribution over . . . , |F|Vl"^'r°}. Note that Wrie) < 

for every e ^ E. Denote w' = w -\- Wr- Let Tyj be an MRCT with respect 
to w and T^/ be that with respect to w' . By Theorem 7, with probability at 
least 1 — 1/(2|1/|), G is strongly min-unique with respect to w' . PARALLEL 
runs FIND-PATH(G,u)', s,t) for an arbitrary pair G F to unambiguously 
test if G is strongly min-unique with respect to w' , and outputs "fail" if it is 
not. PARALLEL then runs APPROX(G, w') to unambiguously output a tree 
T' with 



< + .c^,(r^) (4) 

by invoking Lemma 15 with a sufficiently large constant r such that 8/(9r -|- 
12) < e/2. 

We shall prove that 

c^(T')<(^ + e)-c^(r^), (5) 

which is true by Lemma 16 for sufficiently large \V\. □ 

Sketch of proof of Theorem 20. We omit the simple case where w is identically 
zero and assume without loss of generality that miUgg^; „,(e)-^o ^(^) — 1- Let 
Go be the subgraph of G formed by the zero-weight edges of G. PARALLEL- 
SROCT tests whether Go is a connected spanning subgraph of G in logarith- 
mic space [Rei05] and if so, outputs a spanning tree of Go by calling, say, 
PARALLEL(Go, 0) where denotes the identically zero function. 

Below we assume that Gq is disconnected. The random input to PARALLEL- 
SROCT determines an edge-weight function Wr ■ E ^ Mq" where for each e e E, 
Wr{e) is independently and randomly chosen from the uniform distribution over 
{l/|y|i°,...,|F|V|y|i°}. Let w' = w + Wr. Note that maXeeEWr{e) < 1/|F|^. 
By Theorem 7, G is strongly min-unique with respect to w' with high prob- 
ability PARALLEL-SROCT uses FIND-PATH to determine if G is strongly 
min-unique with respect to w' and outputs "fail" if it is not. Below we as- 
sume that G is strongly min-unique with respect to w' . For each x € V, & 
shortest path tree Tx rooted at x can be output unambiguously by running 
FIND-PATH(G, u)', X, y) for each y £ V. The s.r.c. cost of Tx can be com- 
puted unambiguously by running ROUT-PAIR(T-r, u;', s, t) and accumulating 
the weight of the output path multiplied by r(s) -|- r{t) for all s,t E V. During 
the computation of ROUT-PAIR, each time any bit encoding Tx is needed, it 
is generated on the fly. The spanning tree Tx with the minimum (over x £ V) 
s.r.c. cost with respect to w' is output. The final step in establishing the 
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approximation ratio goes by showing that for every spanning tree T of G, 

v&V 

by the disconnectedness of Go and niingg^ „,(e)^o '"^(c) ^ 1) whereas 



c\^'{T) < m8ix{r{u) + r{v)) 



by Fact 3. The fact that UC C MC C MC^ completes the proof. □ 

Sketch of proof of Theorem 22. We omit the simple case when w is identically 
zero and assume that mingg^; ,„(e)7^o ^('^) — 1- The case where the zero-weight 
edges of G form a connected spanning subgraph of G is dealt with as in the proof 
of Theorem 20, so we may assume that it is not the case. The random input 
to WEIGHTED-2-MRCT determines an edge-weight function Wr : E ^ 
where for each e E E, Wr{e) is independently and randomly chosen from the 
uniform distribution over Let w' = w + Wr- Note that 

maxeeE^^;r(e) < VlV^l"^- WEIGHTED-2-MRCT detects whether G is strongly 
min-unique with respect to w' by running FIND-PATH and outputs "fail" if it is 
not, which occurs with a small probability by Theorem 7. Now, assume that G is 
strongly min-imiquc with respect to w' . The sets , in Theorem 21 where 
dci-) is measured with respect to w' arc computable in unambiguous logarithmic 
space by Theorem 8. For each v E V, let Pv^^^ G SPG(f,si) and Pv^2 ^ 
SPc{v,S2) be the unique shortest paths with respect to w'. By Theorem 8, 
Ti = [j^i-^w' Pv,s\ and T2 = U^g^t"' Pv.s2 unambiguously computable in 
logarithmic space and they are trees by the strong min-uniqueness of G with 
respect to w' . The unique shortest path = (go = si, • • • , Qj,Qj+i, • • • , S2) & 
SFg{si, S2) with respect to w' is also unambiguously computable by running 
FIND-PATH(G, w', si, S2), so is its first vertex qj+i outside of Zf. Theorem 21 
then implies that a tree T satisfying 



c2^(r) < 2c2)(W-2-MRCT^,(G)) 



can be output in unambiguous logarithmic space. The final step in establishing 
the approximation ratio is to show that for every spanning tree T of G, 

c^^\T,si,S2,\)>X+l (6) 

whereas 



&T,Si,S2,X)<)^Tvyu + 



The fact that UC C AfC c ATC^ completes the proof. □ 
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